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Abstract 



Searches for the neutral Higgs bosons h° and A° are used to obtain limits on the Type 
II Two Higgs Doublet Model (2HDM(II)) with no CP-violation in the Higgs sector and no 
additional particles besides the five Higgs bosons. The analysis combines approximately 
170 pb~^ of data collected with the OPAL detector at \fs ~ 189 GeV with previous 
runs at a/s ~ niz and y/s ~ 183 GeV. The searches are sensitive to the h°, A°— >-qq, gg, 
T^T~ and h°— ^>A'^A° decay modes of the Higgs bosons. For the first time, the 2HDM(II) 
parameter space is explored in a detailed scan, and new fiavour independent analyses are 
applied to examine regions in which the neutral Higgs bosons decay predominantly into 
light quarks or gluons. Model-independent limits are also given. 
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1 Introduction 



In this study approximately 170 pb~^ of the dataQ collected by the OPAL detector at 
LEP at 189 GeV centre-of-mass energy are combined with 58 pb~^ of data taken at the 
Z° pole and 53 pb~^ of data at y/s ~ 183 GeV to search for neutral Higgs bosons [1|, 0, Q 
in the framework of the Type II Two Higgs Doublet Model with no CP-violation in the 
Higgs sector and no additional particles besides those arising from the Higgs mechanism 
(2HDM(II)) 0, ^. A model-independent scheme, in which no assumption is made on the 
structure of the Higgs sector, is also analysed. 

In the minimal Standard Model (SM) the Higgs sector comprises only one complex 
Higgs doublet |]l| resulting in one physical neutral Higgs scalar whose mass is a free pa- 
rameter of the theory. However, since there is no experimental evidence for the Higgs 
boson, it is important to study extended models containing more than one physical Higgs 
boson in the spectrum. In particular. Two Higgs Doublet Models (2HDMs) are attractive 
extensions of the SM since they add new phenomena with the fewest new parameters; they 
satisfy the constraints of p ~ 1 and the absence of tree-level flavour changing neutral cur- 
rents, if the Higgs-fermion couplings are appropriately chosen. In the context of 2HDMs 
the Higgs sector comprises five physical Higgs bosons: two neutral CP-even scalars, h° 
and H° (with mh < tmh), one CP-odd scalar, A°, and two charged scalars, H^. 

The most general CP-invariant Higgs potential, having two complex y = 1, SU{2)l 
doublet scalar fields 01 and 02, is given by ^, ^ 

V^(01, 02) = «:i(0l01 - + K2{4<P2 - Vlf + K,[{<P\<P, - vl) + (0^02 - vl)f 

+ «:4[(0l0l)(0^02) - (0l02)(050l)] + ^^,XR<ct>\ct>2) - V.V^f + K^[lm{4>\ct>2)]\ (1) 

where the vacuum expectation values, fj, are non-negative real parameters and the cou- 
plings, real parameters. The physical masses at tree level are given by: 

1 



where 



= 2 [-^11 + -^22 ± V (A^ii - M22? + (2) 
m\ = kq{vI + V2), m\± = K^iyvl + vl), (3) 



Mii = 4(ki + K^)vl + K^vl (4) 

M22 = 4(/t2 + K3)V2 + K^vl (5) 
M12 = {4:K3 + l^5)VlV2. (6) 



The Higgs mixing angle, a, is obtained from 



^ MU-M22 
cos 2a = — , (7j 



sin 2a 



y/{Mu-M22y + ^M 

2.Mi2 

V{Mu-M22y + ^Mj,' 



and the angle /3 is defined as the ratio of the vacuum expectation values, vi and V2, of the 
two scalar fields, tan/3 = V2/V1, with < /5 < 7r/2. 

"'^The searches presented here use subsets of the data sample for which the necessary detector compo- 
nents were fully operational. At ^/s ~ 189 GeV approximately 188 pb~^ were collected and 170 pb~^ 
analysed, varying by ±2% from channel to channel, depending on the detector components required. 



At the centre-of-mass energies accessed by LEP, the h° and A*^ bosons are expected to 
be produced predominantly via two processes: the Higgs-stmhlung process e~ ^h^TP 
and the pair-production process e"'"e~— i>h°A°. The cross-sections for these two processes, 
(Thz and (ThA, are related at tree- level to the SM cross-sections by the following relations 

e+e-->hOzo : a^^ = sm\f3 - a) (9) 
e+e-->hOAO : ahA = cos\(3 - a) X (10) 

where is the Higgs-strahlung cross-section for the SM process e^e^^Hgj^Z*^, and 
^ ~ •^Ah/{-^zh^[12^|/s -|- Azh]} accounts for the suppression of the P-wave cross-section 
near the threshold, with Xij = (1 — m^/s + m^/s)^ — Amfm'j/s'^ being the two-particle 
phase-space factor. 

Within 2HDMs the choice of the couplings between the Higgs bosons and the fermions 
determines the type of the model considered. In the Type II model the first Higgs doublet 
(0i) couples only to down-type fermions and the second Higgs doublet (02) couples only 
to up-type fermions. In the Type I model the quarks and leptons do not couple to the 
first Higgs doublet (0i), but couple to the second Higgs doublet (02)- The Higgs sector in 
the minimal supersymmetric extension of the SM [P, |^ is a Type II 2HDM, in which the 
introduction of supersymmetry adds new particles and constrains the parameter space of 
the model. 

In a 2HDM the production cross-sections and Higgs boson decay branching ratios are 
predicted for a given set of model parameters. The coefficients sin^(/5 — a) and cos^(/5 — a) 
which appear in Eqs. (||) and (|ToD determine the production cross-sections. The decay 
branching ratios to the various final states are also determined by a and (3. In the 
2HDM(II) the tree- level couplings of the h° and A'' bosons to the up- and down-type 
quarks relative to the canonical SM values are 

h°cc:^, h%b:-^, A°cc:cot/3, A°bb: tan/3, (11) 
sm p cos p 

indicating the need for a scan over the range of both angles when considering the different 
production cross-section mechanisms and final state topologies. 

In the analysis described in this paper, detailed scans over broad ranges of these param- 
eters are performed. Each of the scanned points is considered as an independent scenario 
within the 2HDM(II), and results are provided for each point in the (mh, rn^, tan/5, a) 
space. The final-state topologies of the processes (||) and ( p!OD are determined by the 
decays of the Z°, h° and A*^ bosons. Higgs bosons couple to fermions with a strength 
proportional to the fermion mass, favouring the decays into pairs of b-quarks and tau 
leptons at LEP energies. However, with values of a and tan/3 close to zero the decays 
into up-type light quarks and gluons through quark loops become dominant, motivating 
the development of new fiavour independent analyses. 

Section |^ contains a short description of the OPAL detector and the Monte Carlo sim- 
ulations used. The data samples and the final topologies studied are discussed in Section 
1^. The new fiavour independent searches for e"'"e~— *-h°Z° and e"'"e~— >-h°A° are covered 
in Sections ^ and |^, respectively. The model-independent and 2HDM interpretations of 
the searches are presented in Sections || and |^, respectively. In Section ^ the results are 
summarised and conclusions are drawn. 



2 OPAL detector and Monte Carlo samples 



The OPAL detector has nearly complete solid angle coverage and excellent hermetic- 
ity. The innermost detector of the central tracking is a high-resolution silicon microstrip 
vertex detector which lies immediately outside of the beam pipe. Its coverage in polar 
angle0 is | cos6\ < 0.9. The silicon microvertex detector is surrounded by a high preci- 
sion vertex drift chamber, a large volume jet chamber, and z-chambers to measure the 
z coordinates of tracks, all in a uniform 0.435 T axial magnetic field. The lead-glass 
electromagnetic calorimeter and the presampler are located outside the magnet coil. It 
provides, in combination with the forward calorimeter, the gamma catcher, the MIP plug, 
and the silicon-tungsten luminometer |10], geometrical acceptance down to 25 mrad from 
the beam direction. The silicon-tungsten luminometer serves to measure the integrated 
luminosity using small angle Bhabha scattering events The magnet return yoke is 



instrumented with streamer tubes and thin gap chambers for hadron calorimetry and is 
surrounded by several layers of muon chambers. 

Events are reconstructed from charged particle tracks and energy deposits ("clusters") 
in the electromagnetic and hadron calorimeters. The tracks and clusters must pass a set of 
quality requirements similar to those used in previous OPAL Higgs boson searches |jl2| . In 
calculating the total visible energies and momenta, -Evis and P^is, of events and individual 
jets corrections are applied to prevent the double counting of energy of tracks with 
associated clusters 



El- 



A variety of Monte Carlo samples has been generated in order to estimate the detec- 
tion efficiencies for Higgs boson production and background from SM processes. Higgs 
production is modelled with the HZHA generator |]l5l for a wide range of Higgs masses. 
The size of these samples varies from 500 to 10,000 events. The background processes 
are simulated, typically with more than 50 times the statistics of the collected data, by 
the following event generators: PYTHIA {'¥l{l)), grc4f and for the study of the 
systematic errors EXCALIBUR pl (4-fermion processes); BHWIDE (e+e-(7)); KO- 
RALZ H (/i+/i~(7) and r+r-(7)); PHOJET HERWIG |2|, and Vermaseren |3| 
(hadronic and leptonic two-photon processes (77)). The hadronisation process is simu- 
lated with JETSET 
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with parameters described in p4[. The cluster fragmentation 
model in HERWIG is used to study the uncertainties due to quark jet fragmentation. For 
each Monte Carlo sample, the detector response to the generated particles is simulated in 
full detail m. 



3 Data samples and final state topologies studied 

The present study relies on the data collected by OPAL at ^/s ~ mz, 183 and 189 GeV. 
The data collected at the Z° pole provide useful information in 2HDM scenarios where 
the Higgs bosons are light; these data have been extensively analysed in previous OPAL 
publications |]26| , p8| . Higgs search results assuming SM decays from OPAL can be 
found in and |Q for ^ 183 and 189 GeV, respectively. In addition, at y/s ~ 
189 GeV, new flavour independent channels are analysed for the flrst time to explore flnal 
state topologies in which no assumption is made on the quark flavours arising from the 



^ OPAL uses a right-handed coordinate system where the +z direction is along the electron beam and 
where +x points to the centre of the LEP ring. The polar angle, 9, is defined with respect to the +z 
direction and the azimuthal angle, </>, with respect to the horizontal, +x direction. 



Channel h°Z°^ 


Luminosity [pb ^] Data 


Total bkg. 


Efficiency [%] 


= 183 GeV 


bbqq 


54.1 


7 


4.9 ±0.2 ±0.6 


39.2 ± 0.72 ± 1.2 


hhvv 


53.9 





1.56 ±0.13 ±0.18 


47.9 ±0.4 ±0.2 




53.7 


1 


1.3 ±0.1 ±0.2 


41.7±2.5± 1.8 


r+r~qq 


34.2 ±2.1 ± 1.4 


bbe^e^ 


53.7 





0.37 ±0.07 ±0.2 


68.9 ±0.8 ±0.9 


bb/U+/i^ 


53.7 


1 


0.3 ±0.06 ±0.1 


74.6 ±0.7 ±0.7 


^/s = 189 GeV 


bbqq 


172.1 


24 


19.9 ±0.8 ±2.9 


47.0 ± 0.8 ± 1.6 


hhvU 


171.4 


10 


6.9 ±0.5 ±0.6 


42.6 ± 1.1 ± 1.1 


bbr+r^ 


168.7 


3 


4.0 ±0.5 ±0.9 


44.6 ± 1.8 ±2.0 


r+r^qq 


32.8 ± 1.5 ±2.1 


bbe^e^ 


172.1 


3 


2.6 ±0.2 ±0.5 


66.3 ±1.1 ±1.3 


bb/U+/i^ 


169.4 


1 


2.1 ±0.1 ±0.4 


78.3 ± 1.1 ± 1.1 



Table 1: The h^Z" channels: the number of events for the data, the total expected back- 
ground normahsed to the integrated luminosity of the data, and the detection efficiency 
for a Higgs boson decaying only into bb or t^t~ , for typical masses close to the kine- 
matical limits of 85 and 95 GeV at ^/s = 183 and 189 GeV, respectively. Two separate 
efficiencies are shown in the tau channel for the two processes li°— i>bb, TP^t^t^ and 
hP^T^T^ , Z'^^qq. The first error is statistical and the second systematic. 



Higgs boson decays. Detailed descriptions of the flavour independent analyses are given 
in Sections ^ and ^ for the processes e"'"e~^h°Z° and e+e~^h°A°, respectively. 

The channels studied in p9| and using b-tagging, together with those looking 
for r-leptons, provide useful information in regions of the 2HDM(II) parameter space 
where the Higgs bosons are expected to decay predominantly into bb and r+r^ pairs. At 
^/s = 183 and 189 GeV, for the process e"'"e~— >-h°Z° the following flnal states are consid- 
ered: h°Z°— >-bbqq, bbz/i/, bbr+r", r"'"r~qq, bbe"'"e~ and bb/i"*"//". The 2HDM(II) process 
h°Z°— >-A°A°Z°, followed by A°— >bb, is included when kinematically allowed. In addi- 
tion, the 2HDM(II) associated production process, e"'"e~— ^A^h", followed by A°h°— i>bbbb, 
AOhO^bbr+r- (or r+r-bb) and hOAO^AOAOAO^bbbbbb, is studied. 

At ^/s ~ mz the following flnal states are interpreted in the framework of the 2HDM(II): 
h°Z°— >-qqz/z/, qqr+r~, r+r~qq, qqe"'"e~ and qq/i"'"/i", as well as A^h^— ^qqr+r" (or r"'"r~qq), 
and h°A°->A°A°A°— >bbbbbb, if h°^A°A° is kinematically allowed. 

The luminosity, the number of candidate events, the expected SM backgrounds, and 
the efficiencies for each of these h°Z° and h^A*^ channels at 183 and 189 GeV centre-of- 
mass energy are given in Tables |l] and 0, respectively. The signal detection efficiencies 
for the process h'^Z°— i>A°A°Z° can be found in The detection efficiencies quoted 

in Tables |l| and ^ are given as an example for speciflc values of and rriA- When 
scanning the parameter space the efficiency is calculated for each point in the (mh, mA_) 
plane for each of the flnal states considered. In the tau channel two different flnal state 
topologies are studied, in which h°— >bb, Z°— i>r"'"r~ and h"— i>r"'"r~, Z°— i>qq, providing 
different efficiencies, as shown in Table |1|. The dominant contributions to the systematic 
errors on the signal and background efficiencies come from the uncertainty related to 

b-tagging mm. 



Channel A^hO^ 


Luminosity [pb ^] Data 


Total bkg. 


Efficiency [%] 


= 183 GeV 


bbbb 


54.1 


4 


2.92 ±0.2 ±0.5 


50.3 ±0.7 ±2.0 


bbr+r^ 


53.7 


3 


1.50 ±0.1 ±0.2 


44.7 ± 1.6 ± 1.8 


bbbbbb 


54.1 


2 


2.3 ±0.2 ±0.03 


36.0 ± 2.16 ± 1.8 


= 189 GeV 


bbbb 


172.1 


8 


8.0 ±0.5 ±1.4 


48.4 ±0.7 ±3.9 


bbr^r^ 


168.7 


7 


4.9 ±0.6 ±1.6 


45.3 ± 1.5 ±2.3 


bbbbbb 


172.1 


5 


8.7 ± 1.0 ±2.5 


45.4 ±2.2 ±4.3 



Table 2: The h°A° channels: the number of events for the data, the total expected 
background normalised to the integrated luminosity of the data and the signal effi- 
ciency for fmh, mA)=(70 GeV, 70 GeV) and (80 GeV, 80 GeV) in the h°A°-^bbbb and 
bbr"'"r~ channels and for (m^^, mp^)=(60 GeV, 30 GeV) and (70 GeV, 20 GeV) in the 
h^A"— S'A'^A'^A*^— >-bbbbbb channel at y/s= 183 and 189 GeV, respectively. The first error 
is statistical and the second systematic. 

4 Flavour independent searches for e+e~ — ^h'^Z^ 

This section describes the searches for the e"'"e~-^h°Z° process at 189 GeV in 

the following final states: h°Z°— >qqqq and ggqq (the four-jet channel), h°Z°— i>qqz/z/ 
and ggz^z/ (the missing-energy channel), h°Z°— >qqr"'"r~ and ggr+r" (the tau channel), 
h°Z°— >qqe+e~ and gge+e~ as well as h°Z°— >-qq/i+yU~ and gg/x"*"//" (the electron and muon 
channels). A new flavour independent selection has been developed for the four-jet chan- 
nel. The other channels follow closely the analyses described in reference but do not 
make use of b-tagging information. 

4.1 The four— jet channel 

The search in the four-jet channel is a test mass dependent analysis using a binned max- 
imum likelihood method. In order to obtain high sensitivity over a wide range of Higgs 
boson masses, several likelihood analyses are performed. Each of these is dedicated to 
test a specific Higgs mass hypothesis. The test masses are chosen from mh= 60 GeV up 
to mh= 100 GeV in steps of 1 GeV, significantly less than the expected mass resolution 
for the Higgs boson, which is 2 to 3 GeV. Signal Monte Carlo events have been generated 
and reconstructed for each of these test masses. Each likelihood is defined by reference 
histograms made from the background Monte Carlo samples and the signal Monte Carlo 
events generated at the corresponding mass. All data events are then subjected to each of 
the 41 resulting likelihood analyses, and those passing the selection are counted as candi- 
dates for masses in a window of ±0.5 GeV centered on the respective test mass. A single 
event can be a candidate at a variety of different test mass values, and the candidates 
found in the data are not identical for all mass hypotheses between 60 and 100 GeV. 



Correct assignment of particles to jets plays an essential role in separating one of 
the main backgrounds, W"'"W~— >-qqqq, from the signal process, as well as in accurately 
reconstructing the mass of Higgs bosons in signal events. The jet reconstruction method 
is explained in |3^. The initial preselection, designed to retain only events with four 
distinct jets, is unchanged with respect to the four-jet channel in and is independent 
of any mass hypothesis. 

The following selection criteria and the likelihood make explicit use of the mass hy- 
pothesis: 

(1) A kinematic fit which is applied to test the e"'"e~ —>■ h°Z° hypothesis imposes energy 
and momentum conservation. Additionally one of the dijet masses is constrained to 
mz within its natural width and the other dijet mass is constrained to the tested 
mass of the Higgs boson (ZH-Fit). This fit is applied to all six possible jet associ- 
ations and for at least one of these combinations it is required to converge with a 
X^-probability larger than 10~^. The candidate mass is later calculated using the 
jet association which yields the highest x^^probability. 

(2) To discriminate against W'^W" background, the ratio of the matrix element prob- 
ability for Higgs-strahlung |^ and the matrix element probability for W^W~ pro- 



duction as implemented in EXCALIBUR |T8[ is required to be larger than 1.2 x 10 



-4 



When calculating the matrix element probability for Higgs-strahlung it is necessary 
to assign the measured jets to the original partons. Which jets belong to the decay 
products of the Z° and which to the Higgs is determined by the best x^-probability 
of the kinematic fit described in criterion (1). As it remains unknown whether the 
Z° decays into up-type or down-type quarks and which jet belongs to the quark and 
which to the anti-quark, the matrix element is averaged over all these combinations. 
Also, the matrix element probability for W^W~ production is averaged over all pos- 
sible jet-parton assignments. For both matrix element probabilities the four-vectors 
after a four-constraint (4C) fit, imposing energy and momentum conservation, are 
used as input. 

The following six variables are combined with a binned likelihood method [0, with 
one class for the signal and two for the 4-fermion and the 2-fermion backgrounds: 

(a) The logarithm of the ratio of the matrix elements for Higgs-strahlung (MEzh) and 
for W"*"W" production (MEww)- 

(b) The logarithm of the matrix element probability for the Higgs-strahlung process. 
In contrast to the kinematic fits, the matrix element also contains angular informa- 
tion, which allows one to distinguish kinematically between e^e~^h°Z° and Z^Z^^*^ 
production, if mh is in the region of mz- While variable (a) mainly discriminates 
against W"'"W~ events, variable (b) helps to select events compatible with a signal 
hypothesis. The correlations between (a) and (b) are small. 

(c) The logarithm of the x^-probability of the ZH-Fit to the Higgs-strahlung hypoth- 
esis of selection criterion (1). Only the jet association that gives the highest fit 
probability is considered. 

(d) The logarithm of the x^-probability resulting from a kinematic fit (WW-Fit), which 
in addition to energy and momentum conservation forces both dijet masses to be 
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Figure 1: Input variables of the four-jet channel likelihood selection for mh= 80 GeV. 
The OPAL data are indicated by dots with error bars (statistical errors), the 4-ferniion 
background by lighter grey histograms, and the 2-ferniion background by darker grey 
histograms. All MC distributions are normalised to the integrated luminosity of the data. 
The estimated contribution from an 80 GeV SM Higgs boson, scaled up by a factor of 10, 
is shown with a dashed line in each case. 
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Figure 2: Likelihood distribution of the four-jet channel for test masses of (a) 75 GeV 
and (b) 95 GeV. The OPAL data are indicated by dots with error bars (statistical errors), 
the 4-fermion background by lighter grey histograms, and the 2-fermion background by 
darker grey histograms. The contributions expected from a (a) 75 GeV and (b) 95 GeV h*^ 
boson assuming SM cross-section and branching ratios are shown as hatched histograms. 
All MC distributions are normalised to the integrated luminosity of the data. All events 
with a likelihood larger than 0.6 are accepted. 



equal to the mass of the W boson. Only the jet association that gives the highest 
fit probability is considered. 

(e) The difference between the largest and the smallest jet energies in the event. 

(f) The logarithm of an event weight, MEqcdj formed from the tree level matrix 



element for the process e"'"e — >qqqq, qqgg IQ, to reduce QCD background. 



In Figure ^ the distributions of the likelihood input variables are shown for the difficult 
case in which m-w 

The distributions of the final likelihood C^^ are shown in Figure |^ for test masses of 
75 and 95 GeV. Candidates for signal production are required to have C^^ > 0.6 for 
all test masses. In Table |^ the numbers of observed and expected events, together with 
the detection efficiencies for a Higgs signal with a mass of 80 GeV, are given. Figure ^ 
compares the number of candidate events obtained after the likelihood selections for the 
different test masses with the expected background evaluated from Monte Carlo simula- 
tions. In the region of 75 GeV, the number of expected background events rises to more 
than 200. This is due to the presence of W"'"W~ events, where the mass of one of the W 
bosons is constrained to the mass of the Z°, which reduces the other dijet mass by a few 



OPAL 




60 65 70 75 80 85 90 95 100 

m, [ GeV ] 

Figure 3: Number of events that pass the four-jet channel selection at each test mass. 
The number of candidates found in the data is indicated by dots with error bars (sta- 
tistical errors), and the darker and lighter grey histograms correspond to the number of 
events expected from 2-fermion and 4-fermion backgrounds, respectively. The expected 
contribution from a Higgs boson with a mass equal to the test mass, assuming SM cross- 
sections and branching ratios, is shown by the hatched histogram. All MC distributions 
are normalised to the integrated luminosity of the data. Each bin corresponds to a dif- 
ferent analysis of the same data, leading to a strong correlation between neighbouring 
bins. 

GeV compared to the nominal value. The candidate masses are calculated from the 
momenta resulting from a 5C fit requiring energy and momentum conservation and forcing 
one of the dijet masses to mz- Figure | shows the efficiency as a function of the Higgs 
mass for decays to b-quark, c-quark and gluon pairs separately as well as for a mixed 
sample according to SM branching ratios. For mh between 80 and 95 GeV the efficiency 
reaches about 40 to 45% for quarks and 35% for gluons. For small values of mh it drops 
to 25% due to the relatively large amount of initial state radiation that accompanies 
Higgs-strahlung when is considerably lower than the kinematic limit, i.e. i/i— mz. 
For the limit calculation, the efficiencies have been fitted to a polynomial function of mh 
for each flavour separately, and at each mass point the lowest fitted polynomial is used. 
The signal selection efficiencies are affected by the uncertainties given below, expressed in 
relative percentages and shown as an example for a Higgs mass of 75 GeV. The systematic 
errors have been evaluated as follows: for each variable, the distributions obtained from 
the background Monte Carlo samples have been compared to those in the data. Then all 
Monte Carlo events (including the signal samples) have been reweighted so that the mean 
values of data and background distributions become the same and so that the sum of all 
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Figure 4: EfEciency of the four-jet channel selection as a function of mh determined from 
different Monte Carlo samples. The dots represent the efhciency obtained from samples 
with SM branching ratios which also define the reference histograms at each test mass. 
The triangles, squares and crosses correspond to independent samples where the Higgs 
boson decays exclusively to b-quark, c-quark or gluon pairs, respectively. For the limit 
calculation, efficiencies have been fitted to a polynomial function of for each flavour 
separately, and at each mass point the lowest fitted polynomial is used (black line). 

weights is equal to 1. The relative deviations in the number of events which pass the 
selection obtained by reweighting according to the different variables have been added 
in quadrature and amount to 5.5%. The same procedure has been applied to the kine- 
matic likelihood variables, yielding an uncertainty of 2.0%. The uncertainty on the error 
parameterisation of the jet momenta used in the kinematic fits has been evaluated by 
varying the energy and angular resolutions by ±10% , the energy scale by ± 1% and the 
centre-of-mass energy by ±0.3 GeV, each time repeating all kinematic fits. This leads to 
an uncertainty of 6.4%. Since the steps in the test mass are chosen to be smaller than the 
expected mass resolution, the deviation in efficiency due to the interpolation between test 
masses amounts only to 0.7%. The total systematic uncertainty on the signal selection 
efficiency has been calculated by adding the above sources in quadrature yielding 8.7%. 
All of these error contributions have been evaluated for masses of 60 and 90 GeV as well, 
leading to total systematic uncertainties of 9.8% and 7.1%, respectively. The Monte Carlo 
statistical error is about 2%. 

The following uncertainties on the two major background sources are taken into ac- 
count (the first number corresponds to the 4-fermion, the second number to the qq(7) 
background, both for a test mass of 75 GeV): the uncertainty from modelling of the 
preselection cuts is evaluated as described above and amounts to 5.2%/4.5%. For the 
likelihood variables this procedure leads to an uncertainty of 1.2%/1.6%. Varying the 
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Figure 5: Likelihood output of the missing energy channel selection. The OPAL data are 
indicated by dots with error bars (statistical errors), the 4-ferniion background by the 
lighter grey histogram, and the 2-fermion background by the darker grey histogram. The 
estimated contribution from an 80 GeV Higgs with SM cross-section and branching ratios 
is shown as a hatched histogram. All MC distributions are normalised to the integrated 
luminosity of the data. 



error parameterisations, energy scale and centre-of-mass energy for the kinematic fits 
yields an uncertainty of 3.0%/8.1%. Different Monte Carlo generators have been used to 
evaluate the background from 4-fermion events (EXCALIBUR instead of grc4f) and QCD 
events (HERWIG instead of PYTHIA) resuhing in an uncertainty of 2.2%/12.3%. The 
total systematic uncertainty on the residual background estimate amounts to 10.3%. For 
test masses of 60 and 90 GeV the total systematic uncertainties amount to 12.4% and 
8.3%, respectively. The largest Monte Carlo statistical error is 1.3%. 



4.2 The missing— energy channel 

This analysis is nearly identical to a previous one, of which a detailed description can 
be found in with the exception that the b-tagging is not applied. The preselection 



is unchanged. The same kinematic variables used previously, as well as the acollinearity 
angle and the total missing transverse momentum, are combined using a likelihood tech- 
nique. In Figure the resulting signal likelihood distribution is shown for the data, SM 
backgrounds, and an example signal at = 80 GeV. The signal likelihood is required 
to be larger than 0.4 for an event to be selected as a Higgs candidate. 
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Figure 6: Likelihood output of the tau channel selection for events satisfying the 
h°Z°— >-qqr"'"r~ kinematic hypothesis. OPAL data are indicated by dots with error bars 
(statistical errors), 4-ferniion background by the lighter grey histogram, and 2-fermion 
background by the darker grey histogram. The contribution expected from an 80 GeV 
Higgs with SM cross-section and branching ratios is shown as a hatched histogram. All 
MC distributions are normalised to the integrated luminosity of the data. 



The reconstructed mass in selected events is evaluated using a kinematic fit constrain- 
ing the recoil mass to the Z° mass. The numbers of observed and expected events^ are 
given in Table together with the selection efficiencies for an 80 GeV Higgs. The se- 
lection efficiency has been estimated from Monte Carlo samples generated separately for 
b-quark, c-quark and gluon pairs. The efficiency is the lowest and the width of the re- 
constructed Higgs mass is the largest for c-quark pairs and therefore they are used in the 
limit calculation. For an 80 GeV Higgs, the efficiency is (51.0±1.7(stat.)±1.3(syst.))%. 
The efficiency has been optimised for Higgs masses between 70 and 90 GeV. Outside this 
range, the efficiency decreases and reaches about 20% for masses of 60 and 100 GeV. For 
b-quarks, the efficiency is about 5% (relative) higher throughout the whole mass range. 
A total of 47 data events pass the selection, while 44.5±1.4(stat.)±3.0(syst.) events are 
expected from SM background processes. The systematic error has been evaluated as 
in ||2^, but b-tagging related errors have been omitted. 



•^In the calculation of the efficiencies and backgrounds in the missing-energy channel, a 2.5% relative 
reduction has been applied to the Monte Carlo estimates in order to account for accidental vetoes due to 
accelerator related backgrounds in the forward detectors. 



4.3 The tau channel 



The preselection, the tau lepton identification using an artificial neural network, and the 
two-tau likelihood, LrT-, used in this channel are unchanged with respect to the analysis 
described in reference Since b-tagging information is not used, for the final selection 
the likelihood £(qqr"''r^) is used and required to exceed 0.8. Additionally, the 



probability of a kinematic fit, which constrains the invariant mass of the two r's to mz, 
should be larger than 10~^, since this analysis is designed to be sensitive to hadronic Higgs 
boson decays and to Z^^r^r^. In Figure ^the resulting likelihood distributions are shown 
for the data, SM backgrounds, and an example signal at mh = 80 GeV. The numbers of 
observed and expected events are given in Table ^ together with the selection efficiencies 
for an 80 GeV Higgs. The signal detection efficiency has been evaluated for b-quark, c- 
quark and gluon pairs separately and at each mass the lowest value is taken for the limit 
calculation. For an 80 GeV Higgs boson it amounts to (28.7 ± 1.5(stat.)± 2.7(syst.))% 
after all the selection requirements. For lower and higher Higgs masses, the efficiency drops 
to about 20% at mh= 50 GeV and mh= 100 GeV. Two events survive the likelihood cut, 
which can be compared to the expected background of 3.4 ± 0.5(stat.) ± 0.7(syst.). The 
systematic errors quoted above for signal and background are evaluated with the method 
described in [2^, with the contributions from fragmentation and decay multiplicity of 
b-quarks omitted. 



4.4 The electron and muon channels 

The preselection cuts and kinematic likelihood /C are identical to the analysis described 
in ||2^ . Because the present analysis is intended to be independent of the fiavour of the 



Higgs decay products, no b-tagging is applied and /C is used as the final selection variable, 
which should exceed 0.3 for the electron and 0.65 for the muon channel. Figure |^ shows 
the distribution of /C for the electron (a) and the muon (b) channel. 

The signal selection efficiency has been evaluated and fitted for each fiavour separately. 
For the limit calculation, the lowest of these efficiencies at each mass point has been used. 
For an 80 GeV Higgs boson it amounts to (55.4±1.6(stat.)±0.6(syst.))% for the electron 
channel, and (59.3±1.5(stat.)±0.7(syst.))% for the muon channel. In the electron channel, 
the efficiency lies between 50% and 55% for all Higgs masses between 60 and 95 GeV, 
and drops to 27% for a Higgs mass of 100 GeV. In the muon channel, the efficiency is 
between 55% and 70% for all Higgs masses under consideration. 

The numbers of observed and expected events are given in Table ^, together with 
the detection efficiency for an 80 GeV SM Higgs boson. The selection retains 7 events 
in the electron channel and 3 in the muon channel. The total background expectation 
is 4.7±0.2(stat.)±1.4(syst.) events in the electron channel and 4.7±0.1(stat.)±0.9(syst.) 
events in the muon channel. 

The systematic errors quoted above for signal and background are evaluated with 
the method described in ||2^. The largest contribution to the systematic error on the 
background expectation is due to differences between various Monte Carlo generators. 
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Figure 7: (a) The electron and (b) the muon channel likelihood distributions. The OPAL 
data are indicated by dots with error bars (statistical), the 4-fermion background by 
lighter grey histograms, and the 2-fermion background by darker grey histograms. The 
contributions expected from an 80 GeV Higgs with SM cross-section and branching ratios 
are shown as hatched histograms. All MC distributions are normalised to the integrated 
luminosity of the data. 

4.5 Summary of the flavour independent searches for 
e+e-^h^ZO 

The total numbers of candidates accepted after the preselection and the flavour indepen- 
dent likelihood cut, compared to the expected SM backgrounds as well as the detection 
efficiencies for a hadronically decaying Higgs boson with a mass of 80 GeV, are sum- 
marised in Table together with the expected number of signal events in the 2HDM(II) 
for the case of a = 0, tan/3 = 1.0 and mh = 80 GeV. The total number of observed 
candidates for an 80 GeV Higgs boson is 246, while the background expectation amounts 
to 231.3 ± 4.4(stat.) ± 18.0(syst.). 



Cut 


Data 


Total bkg. 


qq(7) 


4-fermion 


Efficiency [%] 


Signal 


Four-jet Channel C = 174.1 pb ^ 


Preselection 


1568 


1521 


379 


1142 


87.4 


31.8 


(1) 
(2) 


696 
648 


649 
600 


117 
116 


532 

484 


65.3 
64.1 


23.8 
23.3 


Likelihood 


187 


174 


28 


146 


35.4 


12.9 


Missing-energy Channel C = 171.8 pb ^ 


Preselection 


111 


101 


18 


83 


63.4 


6.5 


Likelihood 


47 


44.5 


0.6 


43.9 


51.0 


5.3 


Tau Channel C = 168.7 pb ^ 


Preselection 


185 


156 


55 
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49.1 


0.8 


Likelihood 
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3.4 


0.1 


3.3 


28.7 


0.5 


Electron Channel £ = 172.1 pb ^ 


Preselection 


152 


153 


84 


69 


77.8 


1.4 


Likelihood 


7 


4.7 


0.1 


4.5 


55.4 


1.0 


Muon Channel £ = 169.4 pb"^ 


Preselection 


22 


22 


14 


8 


78.8 


1.3 


Likelihood 


3 


4.7 


0.0 


4.7 


59.3 


1.0 



Table 3: The h^Z" channels for the flavour independent analysis: the numbers of events 
selected in the data at -y/s = 189 GeV, expected background from SM processes normalised 
to the data luminosity (given for each channel), the minimal detection efficiencies for a 
Higgs boson with a mass of 80 GeV decaying to quark or gluon pairs and the expected 
number of signal events within the 2HDM(II) for the case of a = 0, tan/? = 1.0 and 
mh = 80 GeV, after the preselection and the flavour independent likelihood cut. In the 
case of the four-jet channel, after preselection, the numbers of events are given for the 
selection with an 80 GeV test mass. 
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Figure 8: Mass distributions for selected data events and SM background expectation 
(a) for all channels and (b) for all channels except the four-jet channel. OPAL data are 
indicated by dots with error bars (statistical) while the background contributions from 
the various channels are shown in different shades of grey. For the four-jet channel the 
90 GeV analysis is shown. 



The mass distributions for candidates found in the data as well as for the background 
expected from the SM for all flavour independent channels are shown in Figure |](a). The 
same mass distributions without the four-jet channel contribution are shown in Figure 



5 Flavour independent search for e+e — A 

We have searched for the process e"'"e~^h°A° in the final states qqqq, ggqq and gggg. The 
signal is characterised by events with four well-separated jets with characteristic invariant 
masses of the dijet pairs originating from the h° and the A°. The dominant background 
is from the process e"'"e~— >-W"'"W~— >qqqq. The second-largest contribution comes from 
e"'"e~— >Z°/7^qq with multiple hard gluon radiation producing a four-jet final state. 

This search is designed to be sensitive over a large portion of the (mh, w^a) plane, 
and the kinematic signatures of signal events depend strongly on mh and niA- For this 
reason, a loose selection is performed first, retaining four-jet hadronic events with partial 
rejection of W"'"W~ and TPTP events. The main discrimination, however, is achieved by 
constraining candidate events to the signal mass hypothesis (mh and wa), and using the 
logarithm of the resulting as the discriminant variable in the limit calculation instead 
of the reconstructed dijet masses. This choice of variable also incorporates naturally the 
measurement uncertainties on the reconstructed masses and simplifies the interpolation 
of its shape as a function of mh and mA. 



5.1 Selection 

Candidate events must first satisfy the requirements of a preselection and then a loose 
selection based on a likelihood variable which is built out of reference distributions of 
reconstructed quantities for events passing the preselection. The following criteria are 
applied ((l)-(4) preselection, (5) selection): 



(1) Each candidate event is required to be classified as a hadronic final state |^ with an 



effective centre-of-mass energy ys' exceeding 150 GeV. The jet resolution parameter 



in the Durham scheme |T3| 1/34 is required to be larger than 0.003 in order to select 



events with four distinct jets. 

(2) Each jet must have at least three tracks. 

(3) The probability of a 4C fit, requiring energy and momentum conservation, must 
be greater than 10"^, to ensure that the mass reconstructions used to isolate the 
signal do not suffer from poor measurement or energy loss from initial state radia- 
tion. 

(4) A 6C kinematic fit is performed requiring energy and momentum conservation and 
also that the invariant masses of the dijet pairs are equal to mw- The 6C fit 
probability of each of the three possible jet combinations is required to be less than 
0.01, to reduce the background from hadronic W"^W~ decays. 



Cut 


Data 


Total bkg. 


qq(7) 


4-fermion 


Efficiency [%] 


Signal 


(1) 


1953 


1885.0 


537.8 


1347.2 


51.2 


14.9 


(2) 


1593 


1532.0 


410.1 


1121.9 


47.9 


13.9 


(3) 


1497 


1446.3 


376.1 


1070.2 


44.3 


12.9 


(4) 


904 


895.0 


329.2 


565.8 


40.0 


11.6 


(5) 


573 


553.2 


238.5 


314.7 


38.5 


11.2 



Table 4: The h^A^ channel for the flavour independent analysis: the numbers of events 
selected in the data at a/s = 189 GeV, expected background from SM processes normalised 
to the data luminosity of 172.1 pb~^, the minimal detection efficiencies for = 30 GeV, 
mA = 60 GeV and the expected number of signal events within the 2HDM(II) for the case 
of a = 0, tan/3 = 1.0, = 30 GeV and = 60 GeV after each step of the selection. 



(5) A likelihood composed of five variables is computed. These variables are the jet 
resolution parameter 7/34, the event-shape variable C obtained from the eigenvalues 
of the sphericity tensor the smallest angle between any two jets in the event, 
the logarithm of the QCD matrix element MEqcd @]' ^"^^ largest proba- 
bility of three 5C kinematic fits constraining energy and momentum and requiring 
the equality of the masses of the two dijet systems (three possible pairings). The 
QCD matrix element used is the maximum of the matrix elements considering all 
possible assignments of observed jets to partons in the e"'"e~— >-qqqq and e"'"e~^qqgg 
processes. The ffist four variables are designed to separate the signal from the qq 
background. The last variable provides rejection of the diboson backgrounds from 
W"'"W~ and TPtP events surviving the 6C fit probability requirement (4). Although 
its use reduces the efficiency for signals with mh = mA, the sensitivity to signals with 
mh 7^ ^A is enhanced. The signal samples used to form reference distributions for 
the likelihood are a mixture of samples in the kinematically accessible region of the 
(mh, mA) plane with mh, mA > 30 GeV. The distribution of this likelihood variable 
for the data, SM backgrounds, and a representative signal with mh = 30 GeV and 
mA = 60 GeV is shown in Figure The likelihood variable is required to exceed 
0.1 for selected events. 

In Table ^ the numbers of events passing the requirements after each step, (1) to 
(5), are given, together with the expected SM backgrounds from 4-fermion and 2-fermion 
processes. The lowest estimated efficiencies, for mh = 30 GeV and mA = 60 GeV, cor- 
responding to the bbbb final state, and the number of expected signal events in the 
2HDM(II) for the case of a = 0, tan/3 = 1.0, mh = 30 GeV and mA = 60 GeV, are shown 
in the last two columns. 

The selection efficiency is estimated with the Monte Carlo simulation at a discrete 
set of reference points in the (mh, wa) plane. The efficiency function is interpolated by 
considering the three closest reference points. A plane in the (mh, mA, efficiency) space 
is formed containing those three points, which allows the efficiency for an arbitrary inter- 
mediate (mh, mA) signal hypothesis to be computed. The interpolated efficiency function 
is shown in Figure for the final-state flavour assignment with the lowest efficiency. 
The efficiency is low near mh = mA = mw because of the difficulty in distinguishing 
hadronic h°A° decays from hadronic W"''W~ decays. The veto of W^W~ events reduces 
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Figure 9: Likelihood distribution for the flavour independent e^e~— i>h°A° channel. The 
OPAL data are indicated by dots with error bars (statistical error), the 4-fermion back- 
ground by the lighter grey histogram, the 2-ferniion background by the darker grey his- 
togram and a representative signal for mh = 30 GeV, m\ = 60 GeV by the dashed 
histogram. All Monte Carlo distributions are normalised to the data luminosity and the 
signal is scaled by a factor of 100. 



the efficiency in the entire (mh, ^a) plane because incorrect jet assignments in h°A° 
events can produce interpretations consistent with the W"'"W~ hypothesis. Conversely, 
the W"*"W~ veto reduces the background of incorrectly-paired W^W~ events everywhere 
in the (mh, mA) plane. The efficiency for low mh and mA is reduced because of the 
requirement that the event should have four distinct jets, which helps to reject the qq 
background. 

After the selection, 573 candidates remain in the data, as compared with the SM 
expectation of 553.2±38.2 events. 



5.2 Discriminant variable 

The invariant masses of jet pairs may be used to separate possible h'^A^ signals from the 
W^W~ and qqgg backgrounds. There are six possible assignments of pairs of jets to h° 
and A° bosons, and there is no constraint from a known mass such as that of the used 
in the e^e~— >h°Z° channels. All six possible interpretations of each selected event are 
tested for consistency with a possible signal. When computing the confidence level for 
excluding a model hypothesis, only a single interpretation of each candidate in the data 
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Figure 10: Lowest efEciency contours for the flavour independent searcli for hadronic 
gyg^fg_ j'jjg efficiency is taken to be zero for mh+mA> 189 GeV. 



and Monte Carlo samples may be used. The choice of jet pairing depends on the mh and 
rriA of the hypothesised signal. 

For each event passing the selection, a 4C kinematic fit constraining energy and mo- 
mentum conservation is performed. For each of the assignments of jet pairs to bosons, the 
reconstructed mh and m\ are computed, along with their covariance matrices. For each 
hypothetical and mp^ considered in the limit computation, each event is assigned the 
jet pairing with the smallest value resulting from the difference of the measured and 
hypothesised mh and mA, and the error matrix of the measurement. The logarithm of 
the smallest is then used as the discriminating variable when computing limits because 
the signal to background ratio depends strongly on the value of logx^- Figure |TI| shows 
the distribution of log for selected data events, the SM expectation, and the signal for 
four mass hypotheses. The signal shown corresponds to the e+e~— *-h°A°— >-gggg process 
because of its poorer mass resolution compared with that obtained for final states with 
quarks. 

For mh = mA = mw, the separation between the signal and the background is poor, 
while for lower values of mh or mA the separation is better. The resolution on the 
reconstructed sum of the dijet masses is approximately 2.4 GeV, while for the difference 
it is approximately 6.2 GeV. The best sensitivity to the signal is in regions of (mh, mA) 
with dijet mass sums different from 2mw- The test mass spacing is determined by the 
model scan grid used when computing the limits - no discretization is introduced within 
the analysis. The scan grid used to compute limits has a finer spacing than the mass 
resolutions on the candidates. All candidates are considered at all test mass hypotheses - 
they simply appear at different locations in the logx^ histogram. The distribution of the 

variable for the signal and backgrounds changes slowly with the test mass hypothesis 



OPAL 




Figure 11: Distribution of logx^ of the mass constraint for the flavour indepen- 
dent e"'"e~— >h*^A° channel for four hypotheses of (m^, mx). The OPAL data are 
indicated by dots with error bars (statistical errors), the 4-fermion backgrounds by 
lighter grey histograms and the 2-fermion background by darker grey histograms. The 
e"'"e~— >h°A°— i>gggg signals are shown with arbitrary normalizations by dashed histograms. 



and is interpolated between Monte Carlo samples generated at different test masses. 

Systematic uncertainties have been considered on the signal and background normali- 
sation and shapes. The e+e~— >-W"'"W~ cross-section is taken to be uncertain at the level of 
2% from a comparison of the predictions of the GENTLE and EXCALIBUR calculations. 
The selection efficiency for e+e~— >W"^W~ events is uncertain at the level of approxi- 
mately 1%, from sensitivity to fragmentation modelling in hadronic W decays and from 
comparisons of the selection variables in data and Monte Carlo [^. The background from 
Z°/7 — >qq(7) has an 11% uncertainty |38|, which includes the uncertainty on the selection 
efficiency and on the four-jet rate in qq events, which is the dominant contribution. The 
4-fermion background from two neutral vector gauge bosons has been estimated using the 
grc4f Monte Carlo generator for the central value, and its uncertainty has been estimated 
by comparing the results obtained with the grc4f and EXCALIBUR generators. Scaling 
these uncertainties by their fractional contributions to the background of this selection and 
adding the results in quadrature yields an uncertainty on the background normalisation 
of 6.9%. Monte Carlo statistics accounts for only a 1% relative error on the background. 



The uncertainty on the signal efficiencies is dominated by the flavour dependence, 
with the highest selection efficiency for the gggg final state. The bbbb and cccc final 
states have very similar selection efficiencies. The lowest signal efficiency at each mass 
hypothesis is used in the limit calculations. 

A more significant effect on the modelling of the signal is the uncertainty in the 
reconstructed mass resolution, as this affects the shape of the logx^ distribution of the 
signal and hence the limits. Similar performances are achieved in Monte Carlo simulations 
of the bbbb and cccc final states, but the gggg final state has on average a positive shift 
of one unit of log relative to the four-quark final states because the resolution is poorer 
for reconstructing masses from gluon jets. The conservative approach of using the logx^ 
distribution of gggg signal final states has been adopted when computing the limits. 



6 Model— independent interpretation 

The results of all the individual search channels at the studied centre-of-mass energies are 
combined statistically to provide 95% confidence level (CL) limits in a model-independent 
interpretation in which no assumption is made on the structure of the Higgs sector. 
The limits are extracted using the same method applied in previous OPAL publications 



Model-independent limits are given for the cross-section of the generic processes 
e"'"e~— > S°Z° and e"'"e^^ S°P°, where S° and P° denote scalar and pseudo-scalar neutral 
bosons, respectively.0 The limits are conveniently expressed in terms of scale factors, 
and 1^91 , which relate the cross-sections of these generic processes to SM cross-sections 
(c.f. Eqs. (i), (0)): 

^sz = (12) 
asP = c'X 4^. (13) 



Figure |T2| shows the 95% CL upper bound for as a function of the S° mass, obtained 
from: 

^sz 



where N^^ is the 95% CL upper limit on the number of possible signal events in the data, e 
is the signal detection efficiency, C is the integrated luminosity, and the sum runs over the 
different centre-of-mass energies of the data and the different channels. In Figure |l^(a) 
only the flavour independent channels described in Section ^ and the channels analysed 
at the pole are used to extract a 95% CL upper limit on s^-BR(S°— >hadrons). In 
Figure 0(b) the SM Higgs branching ratios for the S° are assumed and search channels 
with b-tagging are used. In the region ms < 30 GeV the high energy data (LEP2) 
have little exclusion power while for ms > 50 GeV the Z" data (LEPl) contribute little 
to the determination of the experimental limit. The limit is calculated only for ms 
> 5 GeV, since below this mass value the direct search rapidly loses sensitivity and the 



limit is extracted by a different method ||0|, which makes use of the electroweak precision 
measurements of the Z*^ width and provides an limit of about 0.5 x 10~^. 

The limit on ms for = 1 assuming SM branching ratios is 91 GeV in complete 



agreement with the result obtained by the SM search |Q at \/s~ 189 GeV. A limit of 
75 GeV on ms for = 1 is obtained when assuming a 100% hadronic branching ratio. 

"^Throughout this paper numerical mass hmits are quoted to 1.0 GeV precision. 




Figure 12: Observed and expected upper limits at 95% CL on (a) s^-BR(h^^ hadrons) 
using Havour independent channels, and (b) using all S°Z° search channels with b- 
tagging and assuming SM Higgs branching ratios for the ^ . 



This weaker hmit is partly due to the presence of candidates around ms ~ 80 GeV as can 
be seen from the different behaviour of the observed and expected limit in Figure |l2](a). 

Iso-contours of 95% CL upper limits for in the S° and P° mass plane are shown for 
the processes S°P°— i>qqqq, ggqq and gggg in Figure |13|(a), and for e+e^^ S°P°^bbbb 
and bbr+r^ in Figures |l3|(b) and (c), respectively, assuming a 100% branching ratio into 
the specific final states. The contours are obtained from: 



with N^f being the 95% CL upper limit for the number of signal events in the data. The 
results obtained in Figures |l^(a) and (b) are symmetric with respect to interchanging of S° 
and P*^, while those obtained for r+r^bb are not. For this reason, the results for r+r^bb 
are presented with the mass of the particle decaying into r+r~ along the abscissa and that 
of the particle decaying into bb along the ordinate. The irregularities of the iso-c^ contours 
are due to the presence of candidate events. Along the diagonal, for = 1, a lower bound 
is extracted using the b-tagging channels on the masses at ms = mp ~ mr+r- ~ ""^bb > 78 
GeV at 95% CL. In the hypothesis of S° P° decaying to hadrons with a 100% branching 
ratio, a lower bound of ms = mp > 61 GeV is obtained along the diagonal for = 1. 
Note the small region 30 < mp, ms < 40 GeV in Figures [T^(a) and (b) which is excluded 
by the flavour independent search but not when using only bbbb channels. 
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Figure 13: Upper limits at 95% CL for (a) for the S*^P°^qqqq, ggqq and gggg search 
channel assuming the hadronic branching ratio for both ^ and to be 100%, (b) for 
the S°P°— >bbbb search channel assuming the bb branching ratio for both ^ and P° to 
be 100%, and (c) for the S°P°^bbr+r~ search channel assuming a 100% branching ratio 
for this final state. The invariant masses of the tau-lepton pair and b-jet pair in (c) are 
denoted by m^r sind m^,, respectively. The iso-contour lines are at values of c^ < 1.0, 0.5 
and 0.3 in (a), < 1.0, 0.5, 0.3, 0.15, 0.12 and 0.1 in (b) and < 1.0, 0.5, 0.3, 0.15, 0.1, 
0.008 and 0.003 in (c), respectively. 



7 Two Higgs Doublet Model interpretations 

The interpretation of the searches for the neutral Higgs bosons in the 2HDM(II) is done 
by scanning the parameter space of the modeL Every (mh, rriA, tan/3, a) point deter- 



mines the production cross-section and the branching ratios to different final states. An 
updated version of the HZHA Monte Carlo generator |[T^ that includes the 2HDM(II) 
production cross-sections and branching ratios for Higgs decays has been used to scan 
the parameter space. This generator includes next-to-next-to-leading order QCD correc- 
tions and next-to-leading order electroweak corrections. The branching ratios obtained 



were cross-checked with the results of another generator [M in which QCD corrections 
are computed only up to next-to-leading order. The comparison showed good agreement 
between the results of the two programs. 

The results of all the individual search channels^ at the studied centre-of-mass energies 
are combined statistically to constrain the 2HDM(II) parameter space. Although the 
flavour independent channels supply a unique way to investigate parameter space regions 
where the branching ratio h''— >-bb or A'^^bb is highly suppressed {e.g., low a and tan/3 
regions), they have a poor sensitivity with respect to the b-tagging channels outside these 
regions. The use of b-tagging information substantially reduces the background coming 
from W"'"W~ events and improves the sensitivity to observe Higgs bosons even in regions of 
the 2HDM(II) parameter space where only small branching ratios for h*^— i>bb are expected. 
The expected confidence level is calculated alternatively including only the b-tagged or 
non-b-tagged channels: for each parameter space point, either the flavour independent 
or the b-tagging analysis is then chosen for the extraction of the limits, depending on 
which provides the better expected confidence level. 

The parameter space covered by the present study is: 



• 1 < mh < 100 GeV, in steps of 1 GeV 

• 5 < mA < 100 GeV, in steps of 1 GeV; 
100 < mA < 500 GeV, in steps of 5 GeV; 
0.5 < mA < 2.0 TeV, in steps of 0.5 TeV 

• 0.4 < tan/3 < 58.0, in steps of 1° in (3 

• a = 0, -it/8, -tt/A, -Stt/S and -7r/2 



The values of a are chosen to extend the analysis to the particular cases of maximal and 
minimal mixing in the neutral CP-even sector of the 2HDM(II) {a = —it /A and — 7r/2, 
respectively) and of BR(h°— >bb) = (a = 0). A more complete picture of the model 
is obtained by studying two more intermediate values of a. For tan/3 < 0.4 radiative 
corrections become unstable. Below mA~ 5 GeV the direct search in the e"'"e~— s>h'^A° 
channel cannot be included since the detection efficiency vanishes, and the constraint 
from the total Z° width provides very limited exclusion since the contribution is too 
small. The other two free parameters of the model, mn and mH± , are not scanned in the 
present study. They are fixed at values above the kinematically accessible region at the 
present centre-of-mass energies at LEP2. 

The production of any neutral low mass scalar particle in association with the Z*^ 
was investigated in a previous OPAL publication |^ and, for mh < 9.5 GeV, a mass- 
dependent upper limit on the Higgs boson production cross-section was obtained. This 
limit translates directly into an upper limit on the 2HDM(II) production cross-section for 
mh below 9.5 GeV. Another powerful experimental constraint on extensions of the SM is 

^For the case h^Z"— ^bbr+r" and h'^Z'^^r+r^qq (tau channel) two different efficiencies are appHed, 
according to the final state topologies studied. 



the determination of the total width of the Z° boson at LEP [HDl. Any possible excess 
width obtained when subtracting the predicted SM width from the measured Fz value can 
be used to place upper limits on the cross-section of Z° decays into final states with h° and 



A*^ bosons [Q. An expected increase of the partial width of the Z° is evaluated for each 
scanned parameter space point in the 2HDM(II); if it is found to exceed the experimental 
limit, the point is excluded. The two constraints discussed above are treated together and 
are referred to as Z° width in the rest of the paper, since for low m-^ values most of the 
excluded regions are obtained from the constraints derived from Tz- 

The direct searches for the process e"'"e~— ^-h^'Z^ (e"'"e~— ^-h^A") in the Z° data contribute 
mainly in the m-^< 50 GeV (mh< 60 GeV) region. Since the flavour independent h^'Z^ 
and h°A° analyses have been performed in the mass regions mh> 60 GeV (for the tau and 
missing energy channel, mh> 30 GeV) and m^, mx> 30 GeV, respectively, only b-tagging 
channels using higher energy data are applied below these masses; however these channels 
have no detection efficiency for mh< 30 GeV. The flavour independent analyses provide 
exclusion for the whole tan/3 range and for the tan/?<l regions for a = and a = — vr/8, 
respectively. In Figures 0(a-e) the excluded regions in the (mh, mp^) plane are shown 
for the flve chosen values of a, together with the calculated expected exclusion limits. A 
particular (mh, niA, a) point is excluded at 95% CL if it is excluded for all scanned values 
of tan/3. Different domains of tan/3 are studied and described below: a) 0.4 < tan/3< 58.0 
and b) 0.4 < tan/3 < 1.0 or 1.0 < tan/3< 58.0, for which enlarged excluded regions are 
obtained. 

a) 0.4 < tan/3< 58.0 (darker grey area): 

• The poor sensitivity of the channels below mh < 10 GeV causes a sharp cut in the 
exclusion plots at this mh for a > —7i/2. The exclusion in this region is extracted 
from the total width of the boson, as explained above. Both the h°Z° and h^A^ 
production processes contribute to the natural width of the Z°. While an excess, 
induced by the h°Z° process, extends the exclusion region to any value of ttia, the 
exclusion provided by the h'^A'' process is kinematically limited to the region where 
mA + mh < TUz- The contribution of the h°Z° production cross-section to the Z° 
width depends on the argument (/3 — a), and it becomes large enough for this process 
alone to provide exclusion in different tan/3 domains for the a values considered. 

• For a = and —vr/8, most of the exclusion is provided by the channels at y/s = mz, 
where no b-tagging was applied. In fact, the flavour independent analyses at ^/s ~ 
189 GeV have a limited sensitivity because of the presence of the W"'"W~ background 
events. The line at mh ~ 57 GeV in Figure |14|(b) is a result of the Z° data kinematic 
constraint. 

• The presence of candidates in the four-jet and r+r^ b-tagging h^A'' channels at 
y/s = 189 GeV at mh, mA ~ 80 GeV, due to the W^W~ background, is clearly 
reflected in Figure |l5(c). The unexcluded region at (mh, mA) ~ (60, 90) GeV in 
0(d) is due to the presence of candidates in the four-jet b-tagging h°A° channel. 

• The shape of the exclusion plot in Figure |l^(e) for mh < 35 GeV is related to the 
kinematical constraint on the h°A° production in the Z° data, which for a = —7c/2 
and large tan/3 is the only allowed process, since the h°Z° production cross-section 
vanishes when /3 — a ^ vr. For mh > 35 GeV, the high energy data open a new 
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Figure 14: Excluded regions in the fmh, mj>^) plane, (a)-(e), for a = 0, —tc/S, — 7r/4, 
— 37r/8 and —it/2, respectively, together with the expected exclusion limits. A particular 
fmh, rriA, Oi) point is excluded at 95% CL if it is excluded for all scanned values of tan/?. 
Three different domains of tan/? are shown: the darker grey region is excluded for all 
values 0.4 < tan/5< 58.0; additional enlarged excluded regions are obtained constraining 
0.4 < tan/3 < 1.0 (lighter grey area) or 1.0 < tan/3< 58.0 (hatched area). Expected 
exclusion limits are shown for 0.4 < tan/?< 58.0 (dashed line). 



kinematic region and are able to exclude large (mh, mA) areas, as can be seen by 
the sharp line in Figure p^(e) . 



• The (mh, points below the semi- diagonal defined by mh> 2mA, for which 
the process h"— >A°A° is kinematically allowed, can only be excluded for tan/3 > 
0.5 values by the high energy channels. In fact, for very low values of tan/? the 
branching ratio for A°— s>bb vanishes, causing unexcluded regions in Figures |T^(c), 
0(d) and |l^(e), which are excluded by the Z° data fiavour independent analyses 
below mh ~ 60 GeV. 

b) 0.4 < tan/9 < 1.0 (lighter grey area) and 1.0 < tan/?< 58.0 (hatched area): 

• As discussed above, as a consequence of the variation of the h^Z" production cross- 
section with tan/5 in the mh < 10 GeV region, for a = — 7r/8 in Figure |l^(b), the 
mh > 7 GeV region is excluded for all values of mA in the tan/5 < 1.0 domain. For 
a = — 7r/4 and — S/Svr, the mh < 10 GeV region is excluded for all values of mA 
only in the tan/5 < 1.0 domain. 

• At a = and a = — vr/S and small values of tan/5 the production cross-section 
for the process e"'"e~— *h°Z° is highly suppressed. For mh > 40 GeV, constraining 
tan/? > 1.0, larger excluded regions are obtained, as can be seen in Figures |T^(a) 
and (b) (hatched areas). 

• The presence of candidates in the four-jet and t~^t~ b-tagging h^A^ channels at 
\/s= 189 GeV at mh, mA ~ 80 GeV, corresponding to the W+W" background, 
is clearly reflected in Figures |l5(c), (d) and (e). For a = —n/A and a = —Stc/S 
this region is unexcluded even for tan/? > 1.0, while for tan/? < 1.0 it is excluded 
due to a large expected production cross-section. For a = —n/2, the production 
cross-section becomes small and this domain is unexcluded for tan/5 < 1.0, as can 
be seen in Figure 0(e). 



In Figure |T5| the excluded regions in the (mh, mA) plane independent of a are given 
together with the calculated expected exclusion limits. A particular (mh, mA) point is 
excluded at 95% CL if it is excluded for all scanned values of tan/5 and a. Different 
domains of tan/? are shown: 0.4 < tan/? < 58.0 (darker grey area), 0.4 < tan/? < 1.0 
(lighter grey area) and 1.0 < tan/5 < 58.0 (hatched area), for which enlarged excluded 
regions are obtained. The rectangular region 1 < mh < 44 GeV for 12 < mA ^ 56 GeV 
is excluded at 95% CL independent of a and tan/?. The cross hatched region shows the 
exclusion provided by the constraints on the width of the Z" common to all the scanned 
values of a and tan/5. 

In Figures p!6|(a-e) the excluded regions in the (tan/?, mh) plane are shown for the five 
chosen values of a, together with the calculated expected exclusion limits. A particular 
(mh, tan/3, a) point is excluded at 95% CL if it is excluded for all scanned values of mA. 
There are two regions shown, the whole domain 5 GeV < mA< 2 TeV (darker grey area) 
and a restricted domain for which 5 < mA < 60 GeV (lighter grey area). The exclusion 
contours for mA < 60 GeV are larger for all a values, and entirely contain the 5 GeV < 
"^A< 2 TeV excluded areas. 

In Figures f6Ka) and (b), the unexcluded regions at low tan/? and mh < 10 GeV refiect 
the behaviour in Figures |T5(a) and (b) in the (mh, mA) plane for the same values of a, 
for mh < 10 GeV and mA > 60 GeV. Similarly, the remaining three values of a show a 
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Figure 15: Excluded (mx, m^) region independent of a, together with the expected 
exclusion limit. A particular (m^, m^) point is excluded at 95% CL if it is excluded for 
0.4 < tan/? < 58.0 (darker grey region), 0.4 < tan/? < 1.0 (lighter grey region) and 1.0 
< tan/? < 58.0 (hatched region) for any a. The cross-hatched region is excluded using 
constraints from Fz only. Expected exclusion limits are shown as a dashed line. 



complementary behaviour at large tan/?, reflecting the corresponding exclusion regions in 
the (mh, mp^) plane, Figures |T^(c), (d) and (e). In the mA < 60 GeV contours the same 
regions are obviously excluded. The h°Z° production cross-section contribution to the Z° 
width is not large enough to exclude the small region tan/?> 10, 7 < m-h < 10 GeV in 
Figure |16|(b). In Figure |l^(e) the two unexcluded regions in the lighter grey area (5 < 
mA < 60 GeV) at mh ~ 10 and 35 GeV are projections of the two regions unexcluded at 
the same mh values in Figure ^^e) (darker grey area), for mA < 10 GeV and mA < 60 
GeV, respectively. 

In Figure |1^ the excluded regions in the (mA, tan/3) plane are shown for different 
values of a, together with the expected exclusion limits. A particular (mA, tan/?, a) point 
is excluded at 95% CL if it is excluded for all scanned values of mh. There are three 
regions shown, corresponding to different mh domains that are subsets of one another, 
namely: 1 < mh < 90 GeV (darker grey area), 1 < mh < 75 GeV (lighter grey area) and 
1 < mh < 60 GeV (hatched area). The lower the mh upper value analysed, the larger the 
excluded (mA, tan/5) region. The unexcluded regions for mA > 60 GeV in Figures p!7|(b- 
e) gradually increase with decreasing a. There is an exact correspondence between the 
largest excluded tan/3 value and the excluded regions at mh < 10 GeV in the (mh, tan/5) 
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Figure 16: Excluded regions in the (tan/?, 171^) plane, (a)-(e), for a= 0, — vr/S, — 7r/4, 
—3n/8 and —7i/2, respectively, together with the expected exclusion limits. A particular 
(rriY,, tan/3, a) point is excluded at 95% CL if it is excluded for all scanned values of mp^. 
The two regions shown correspond to the whole domain 5 GeV < mx< 2 TeV (darker 
grey area) and a restricted domain for which 5 < < 60 GeV (lighter grey area). The 
exclusion regions for < 60 GeV entirely contain the 5 GeV < mA< 2 TeV excluded 
areas. Expected exclusion limits are shown for 5 GeV < mA< 2 TeV (dashed line). 
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Figure 17: Excluded regions in the (mx, tan/?j plane, (a)-(e), for a= 0, — vr/S, — 7r/4, 
— Svr/S and —-njl respectively, together with the calculated expected exclusion limits. A 
particular (m^, tan/3, a) point is excluded at 95% CL if it is excluded for all scanned 
values of nii^. The three contours correspond to 1 < mh < 90 GeV (darker grey area), 
1 < rriYi < 75 GeV (lighter grey area) and 1 < rrih < 60 GeV (hatched area). Expected 
exclusion limits are shown for 1 < mh < 90 GeV (dashed line). 



projections in Figures p!B|(b-e) . The line at mA ~ 30 GeV for tan/3 > 1.0 in the exclusion 
region for mh < 90 GeV in Figures 0(d) and (e) corresponds to the horizontal line at 
mA = 30 GeV in the tan/? > 1.0 contour in Figures |l5(d) and (e). The small island at 
35 < rriA < 39 GeV and 1.0 < tan/? < 2.0 in Figure ^]\{e) is the reflection of a candidate 
in the h^A*^— s>r+r^bb channel at ^/s= 183 GeV with mA ~ 37 GeV and mh ~ 80 GeV. 

8 Conclusions 

A general analysis of the 2HDM(II) with no CP-violation and no extra particles besides 
those of the SM and the five Higgs bosons has been performed for the first time. Large 
areas of the parameter space of the model have been scanned. In the scanning procedure 
the dependence of the production cross-sections and branching ratios on the angles a and 
/?, calculated with next-to-next-to-leading order QCD corrections and next-to-leading 
order electroweak corrections, has been considered. 

In addition to the standard OPAL b-tagging analyses, new flavour independent chan- 
nels for both the Higgs-strahlung process, e"'"e~— s>h°Z°, and the pair-production process, 
e"'"e~— >-h°A°, have been analysed, providing access to regions of parameter space in the 
2HDM(II) where h" and A" are expected to decay predominantly into up-type light quarks 
and gluons {e.g. a ~ 0). 

OPAL data collected at ^/s ~ mz, 183 and 189 GeV have been interpreted both in the 
context of the 2HDM(II) and in a model-independent approach where both SM branching 
ratios and 100% hadronic branching ratios were assumed. 

The 2HDM(II) parameter space scan, for 1 < mh < 100 GeV, 5 GeV < mA < 2 TeV, 
—tt/2 < a < and 0.4 < tan/3 < 58.0, leads to large regions excluded at the 95% CL in 
the (mh, mA) plane as well as in the (mh, tan/3) and (mA, tan/3) projections. The region 
1 ^ "^h ^ 44 GeV and 12 < mA < 56 GeV is excluded at 95% CL independent of a and 
tan/3 within the scanned parameter space. 

In the model-independent approach for e+e^— S°Z°, lower bounds at 95% CL are 
obtained for = 1 of 91 GeV assuming SM branching ratios and 75 GeV assuming 
100% hadronic branching ratios. In the case of the generic processes S'^P*^ bbbb and 
gopo _^ T^T~ , a lower bound at 95% CL of mg = mp > 78 GeV is extracted along the 
diagonal for = 1 assuming 100% branching ratios for the individual final states, while 
assuming the S° P" hadronic branching ratios to be 100% gives ms = mp > 61 GeV. 
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